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ABSTRACT 

We observed M51 at three frequencies, 1.4 GHz (20cm), 4.9 GHz (6cm) and 
8.4 GHz (3.6cm), with the VLA and the Effelsberg 100m telescope to obtain 
the highest quality radio continuum images of a nearby spiral galaxy. These 
radio data were combined with deconvolved Spitzer IRAC 8 jjim and MIPS 24 //m 
images to search for and investigate local changes in the radio-IR correlation. 
Utilizing wavelet decomposition, we compare the distribution of the radio and IR 
emission on spatial scales between 200 pc and 30kpc. We show that the radio- 
IR correlation is not uniform across the galactic disk. It presents a complex 
behavior with local extrema corresponding to various galactic structures, such as 
complexes of H II regions, spiral arms and interarm filaments, indicating that the 
contribution of the thermal and non-thermal radio emission is a strong function 
of environment. In particular, the relation of the 24 //m and 20 cm emission 
presents a linear relation within the spiral arms and globally over the galaxy, 
while it deviates from linearity in the interarm and outer regions as well in the 
inner region, with two different behaviours: it is sublinear in the interarm and 
outer region and over-linear in the central 3.5 kpc. Our analysis suggests that the 
changes in the radio/IR correlation reflect variations of ISM properties between 
spiral arms and interarm region. The good correlation in the spiral arms implies 
that 24 yum and 20 cm are tracing recent star formation, while a change in the 
dust opacity, 'Cirrus' contribution to the IR emission and/or the relation between 
the magnetic field strength and the gas density can explain the different relations 
found in the interarm, outer and inner regions. 
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Subject headings: galaxies: individual (M51a, NGC5194) - galaxies: ISM - infra- 
red: galaxies - radio continuum: galaxies 



Introduction 



The observed empirical correlation between the total radio and far infrared (FIR) emis- 



sion of galaxies is very tight (e.g. Helou et al. 1985; de Jong et al. 1985). This relation 



holds over five orders of magnitude in luminosity, from nearby galaxies (normal spirals, 



starburst galaxies, including dwarf, irregular and Seyfert galaxies; Condon 1992 Yun et al. 



2001t |Bell||2003[ ) to distant objects at z>3-4 ( [Appleton et al.||2004[ [Seymour et al.||2009 
Sargent et al. 2010a|[b ). This correlation is interpreted as the result of recent massive star 
formation. Ionized UV photons emerging from young massive stars (>8Mq) heat the sur- 
rounding dust, which then reradiates at IR wavelength. These stars evolve to supernovae 
(SN), whose remnants accelerate cosmic ray (CR) electrons which then emit non-thermal 
(NT) radio synchrotron emission while propagating through the interstellar magnetic field 
and losing energy. The radio continuum emission in normal star-forming galaxies at 20 cm 
wavelength consists of the dominant NT emission, attributed to synchrotron emission from 
cosmic ray electrons, and the less important thermal free- free emission emerging directly from 



Hll regions (Condon 1992). At 3.6cm, the ratio of thermal fraction to the radio continuum 



emission can be up to 50% (Tabatabaei et al. 2007b). 



The IR radiation is primarily emitted from dust heated by H ll regions and cool dust 
excited mainly by the interstellar radiation field. While the thermal radio continuum and 
warm dust emission can be both directly associated with Hll regions, the IR and NT radio 
emission mechanisms involve very different physical processes and time scales. Thus the 
origin and physics of the NT radio-IR correlation remain still under debate. Even though 
this correlation is still poorly understood, it is used to calibrate the radio emission as a 



tracer of recent star formation in local and high redshift galaxies (Condon 1992 Yun & 



CariUi 


2002 


Bell 


2003 



hence thought to be a good estimator for SF in galaxies, as radio interferometers offer much 
higher spatial resolution than space-based IR telescopes. The local radio/FIR correlation 
within galaxies has been investigated ( Beck fc Golla|[l988 Hoernes et al.||1998 Hughes et al. 



2006a Murphy et al. 2006 Tabatabaei et al. 2007a) in order to probe in detail the origin of 



this correlation and the mechanisms at work. In particular the determination of the spatial 
scales on which the radio-IR correlation breaks down can provide strong constraints on the 



different models explaining this correlation (Volk 1989 Helou & Bicay 1993 Niklas & Beck 



1997). These studies showed that while the peaks of the radio and IR emissions are spatially 



3 



close (Hoernes et al. 1998), the radio map seems to be a smeared version of the IR map 



( Murphy et aL||2006 ), reflecting the CR electrons' diffusion in the interstellar medium (Bicay 



& Helou 1990). 



In this paper we investigate the radio continuum emission and its correlation with the 
IR emission within M51 at high spatial resolution, down to ~200pc. Throughout this 
paper, M51 refers to the interacting galaxy pair while M51a and M51b refer, respectively, 
to NGC5194 and the companion NGC5195. The aim of this study is to probe on what 
spatial scales the radio-IR correlation holds and the role played by the environment for this 
relation (e.g. arm/interarm regions). In particular, an important point to constrain is how 
well the radio emission is correlated with SF regions, keeping in mind the ultimate goal of 
constraining the NT radio continuum as tracer of recent star formation events. M51a is an 
ideal target to conduct such a study. It is nearby (8.2 Mpc, which corresponds to ~40pc per 
arcsec, TuUy 1988, with Hq = 70 kms~^ Mpc~^), has well-defined spiral arms and is almost 
face-on (disk inclination of ~20°). The high disk surface brightness allows one to study in 
detail arm/interarm variations. It is also one of the best studied galaxies at almost every 
wavelengths. 

We obtained VLA data at three wavelengths: 20 cm (1.4 GHz), 6 cm (4.9 GHz) and 
3.6 cm (8.4 GHz) at high resolution (~2") and high sensitivity (rms < 25/iJy/beam) to 
probe the distribution of the radio continuum emission across the disk of the galaxy. We 
compare the radio continuum maps with SPITZER MIPS 24 /xm and IRAC 8 /im images 
and examine the radio-IR correlation down to 200 pc, using the wavelet transformation and 



cross-correlation method introduced by Frick et al. (2001). This analysis technique has been 



successively applied to NGC6946 by Frick et al. (2001), M33 by Tabatabaei et al. (2007a) 



and the Large Magellanic Cloud by Hughes et al. (2006a) and allows us to investigate the 



radio and IR emission and their correlation as a function of the different galactic structures. 

The paper is organized as follows: In Sect. [2] we present our radio observation and their 
data reduction plus the ancillary data used here. We describe the radio continuum emission 
of M51 in Sect. |3} In Sect. |4]we conduct the structural analysis of the radio and IR images 
of M51a using continuous 2D wavelet transformation and we investigate the results of this 
analysis in the context of the radio-IR correlation in Sect. |5| 



2. Observations and Data reduction 



In the following we describe the datasets used in this paper. An overview is provided 
in Table [H 
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2.1. Radio observations 

M51 was observed with the Very Large Array (VLA) at 20 cm (1.4 GHz) in the A, B, 
C and D configuration, at 6 cm (4.9 GHz) in the B, C and D configuration and at 3.6 cm 
(8.4 GHz) in the C and D configuration. At 3.6 cm, the primary beam FWHM is 5'A, thus 
two pointings were necessary to cover the full disk of M51 with a size of 10'. The VLA 
data were taken between 1988 and 2005 as part of several projects, corresponding to a total 
integration time of ~170hrs on-source spread over all three frequencies. Table |2] provides a 
summary of the different observation runs. For the two most compact configurations (C and 
D) the continuum correlator mode was used, while observations in A and B configuration 
were done in multichannel continuum mode, using the four-IF spectral line mode. Both 
setups use two intermediate frequencies (IF). The central frequencies and corresponding 
bandwidths are listed in Table |2j During all observations, the quasar 3C286 was observed 
for fiux calibration. For the 20 cm observations in C configuration, the quasars 3C48, 3C138 
and 3C138 were also available as additional fiux calibrators. The different projects used 
the following sources as phase calibrators: 1335+457, 1349+536, 1409+524, 1219+484 and 
137+434. Depending on the frequency and configuration, the phase calibrator was observed 
every lOmin to 45min. 

Data reduction was done using standard AIPS routines. Each day of observation has 
been calibrated separately. We performed standard flux calibration and atmospheric correc- 
tion for the uv-data of each observation day. Before and after those calibrations, the data 
were inspected for RFI (for each IF and polarization) and flagged in consequence. For the 
data observed in the multichannel correlator mode (20 cm in A and B conflguration and 6 cm 
in B conflguration), calibration and flagging were done on the pseudo-continuum channel and 
then applied to the 'Line' data. Bandpass calibration was also performed on the 'Line' data 
using the flux calibrator and the pseudo-continuum channel. Finally we extract the source 
from the calibrated data with the task SPLIT which produces single-source flies and applies 
the calibration. All days of observation in one conflguration were then combined together 
with the task DBCON. At this point, we have one fully calibrated single-source data set for 
each conflguration and each frequency. 

For each frequency, we combined the fully calibrated uv data in the different conflgu- 
rations. First, we converted the coordinates of the oldest uv data from the B1950 to J2000 
coordinates system with the tasks REGRD and UVFIX (20 cm data in C and D conflguration 
and 6cm data in D conflguration). Then, the intermediate frequencies being different for 
the different projects, the IFs were splitted with the task SPLIT for each conflguration, and 
then recombined with the task DBCON into a single-channel data set. In the case of uv data 
obtained in the multichannel continuum correlator mode (conflgurations A and B at 20 cm 
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and 6 cm) , this step was done after all channels were splitted and recombined into a single 
data set. Therefore, at each frequency, a single data set was obtained for all configurations 
which were then combined in chronological order (the older observations first) with DBCON 
to obtain a final uv data set at a given reference frequency. 

In a next step, the radio maps are produced using the task IMAGR. As the 20 cm and 
6 cm data cover a large field-of-view, we first divided the field into 13 and 55 facets, of 
4096x4096 pixels and 2048x2048 pixels at 6 cm and 20 cm, respectively. Each facet at 20 cm 
and 6 cm (respectively, each pointing at 3.6 cm) was then imaged, with the multi-resolution 



CLEAN in AIPS (Rich et al. 2008 Greisen et al. 2009), using several circular Gaussian 



model sources, in order to optimally deconvolve and image the extended faint emission. At 
6 cm and 3.6 cm, we choose a number of circular Gaussian models equal to the number of 
observed configurations: 2 Gaussian models at 3.6 cm, of width equal to 0" (point source) 
and 6'.'3 (with a flux density cut-off of 35 and 40mJy/beam, respectively), 3 Gaussian models 
at 6 cm, of width 0", 6'.'3 and 20" with flux density cut-offs of 0.05, 0.06 and 0.2mJy/beam, 
respectively. At 20 cm, 4 Gaussian models of width 0", 4'.'7, 15" and 47" were chosen (with a 
respective flux density cut-off of 0.03, 0.3, 3 and 30 mJy/beam). The widths of the Gaussians 
are multiple of the robust beam size (respectively 2'.'4, 2"0 and 1"5 at 3.6cm, 6cm and 20cm 
respectively), separated by a factor of a/IO, thus giving a factor of 10 between the beam 
area of successive Gaussians. For the data at 20 cm and 6 cm we used respectively a robust 
weighting of ROBUST = 1 and ROBUST = 0, and at 3.6 cm we deconvolved with uniform 
weighting for maximal resolution (ROBUST= —4). After the deconvolution, we combined 
the facets/pointings using the task FLATN to construct the flnal images at 20 cm, 6 cm 
and 3.6 cm. This task performs also a correction for the primary beam attenuation while 
combining the two pointings at 3.6 cm. During this correction, pixels outside the 10% power 
radius of the primary beam were blanked. 

At 6 cm and 3.6 cm, the largest scale structures detectable (5' and 3', respectively) are 
smaller than the size of the galaxy system (10'). This implies that at these wavelengths 
the VLA observations can miss extended emission from M51. To correct for missing short- 



spacings, we use the same method as Tabatabaei et al. (2007b), which consists of combining 



the flnal VLA CLEANed images with single dish data obtained with the Effelsberg 100m 



telescope (Fletcher et al. 2010). First, the VLA maps were convolved to a circular beam, 
2" at 6 cm and 2.4" at 3.6 cm, the corresponding Effelsberg beam size being 180" and 90", 
respectively. We subtracted strong point sources in the VLA images to avoid bias in the 
combined total flux and the VLA and Effelsberg maps were combined in the uv-plane with 
the task IMERG in AIPS. The only free parameter of this task is the choice of the uv range 
of the interferometric data that should be flUed by the single-dish values. We choose the 
uv range such that the integrated flux densities of the resulting merged map and of the 
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Effelsberg map are equal. The overlap uv-range is O.SlkA to 0.71kA at 6 cm and l.OlkA 
to 1.33kA at 3.6 cm. The VLA point sources were then added back to the short-spacing 
corrected (SSC) maps. The total flux densities of the resulted SSC and the corresponding 
single dish maps are in agreement withing 0.1%. Short-spacing correction is not necessary at 
20 cm, since observation from the compact D configuration of the VLA allows one to detect 
structures with sizes up of 15', larger than M51. 

Finally, we correct the 20 cm and SSC 6 cm maps for primary beam attenuation with the 
task PBCOR, with the same cut-off value of 10% of the primary beam radius used for the 
3.6 cm image. The final 3.6 cm SSC, 6 cm SSC and 20 cm maps corrected from primary beam 
attenuation have 4096x4096 pixels with a pixel size of 0'.'4, 0'.'2 and 0'.'3 at 3.6 cm, 6 cm and 
20cm, respectively, and a final resolution of 2'.'4x2'.'4 at 3.6 cm, 2"x2" at 6 cm and r.'4xr.'3 
at 20 cm. The r.m.s. noise level is similar for the 20 cm (cr„o«se=ll /iJy/beam) and 6 cm data 
(o"nojse=16 /iJy/beam), while it is higher in the case of the 3.6 cm map ((T„ojse=25 /iJy/beam). 



2.2. Spitzer infrared data 



For comparison with our high resolution radio maps, we used the Spitzer MIPS 24 /xm 
map from the fifth Spitzer Infrared Nearby Galaxies Survey (SINGS) data deliver}|^ and the 
8 /im PAH emission map, which corresponds to a stellar continuum subtracted 8 /im image. 



created by Regan et al. (2006). The description of the observations and reduction can be 



found in Kennicutt et al. (2003) and Calzetti et al. (2005). Briefly, the data reduction of these 



images used the standard MIPS pipeline ( Gordon et al.|[2005 ). In addition, some extra steps 
were performed to improve the image quality. For the MIPS 24 /im image, fiat-fielding was 
first performed, hot and dead pixels were removed and latent images from bright sources 
and the background were subtracted. Cosmic rays were removed while the mosaic was 
created. For the IRAC images additional steps include correction from geometric distortion 
and rotation, image offsets and bias drift, cosmic ray removal and background subtraction. 



Finally the stellar continuum was subtracted (Regan et al. 2006). The final resolution at 



24 /xm is 5'.7 (Fig. |l| top left) and it is 2" for the 8 /im PAH emission map (Fig. [T| bottom 
left). 

As the MIPS point spread function (PSF) has strong secondary lobes and Airy rings 
which would affect a direct comparison to the radio images, we applied the HiRes deconvo- 
lution algorithm ( Backus et al.|[2005 ) to the IR maps. In addition, this deconvolution allows 



for an increase of the resolution of up to a factor of two in FWHM (Velusamy et al. 2008). 



ihttp://irsa.ipac.caltech.edu/data/SPITZER/SINGS/ 
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Therefore we can analyse our radio and IR data and investigate the correlation between the 
two without degrading our radio data by convolution with the MIPS PSF as it has been 



done in a previous study (Murphy et al. 2006). HiRes is based on the maximum correlation 



method (MCM), an extension of the Lucy- Richardson deconvolution algorithm and has been 



well tested on IRAC as well as MIPS data (Backus et al. 2005 Velusamy et al. 2008). The 



deconvolution has been carried out with 50 iterations on the fully reduced SINGS MIPS mo- 
saic at 24 nm and on the PAH IRAC 8 fim maps, using the point response functions (PRF) 
provided by the Spitzer Science Center (SSC). The resulting HiRes deconvolved images are 
shown in Fig. [l| along with the initial SINGS images at 24 yum (top left and top right, re- 
spectively) and 8 fim (bottom left and bottom right, respectively). The HiRes deconvolution 
provides a FWHM in the range of 0'.'9 - 1'.'5 and 2'.'2 -2'.'6 for the 8 fim and 24 /im images, 
respectively, as obtained by fitting Gaussian models to several point sources in the vicinity 
of M51. The spread in the FWHM is caused by a lack of very bright stars in the FOV. We 
adopt a final resolution of I'.'lx 1'.'3 at 8 fim and 2'.'3x2'.'4 at 24 /im. 

The benefit of applying the HiRes deconvolution is improved visualization of spatial 
morphology by enhancing resolution and removing the contaminating sidelobes. Removing 
the diffraction lobes in the Spitzer images is more critical than resolution enhancement for 
our analysis of the radio/IR correlation. The sharper inner edges of the spiral arms become 
clearly visible in the HiRes image. The 24 /xm HiRes map reveals in detail the filamentary 
structure in the inter-arm regions and allows us to separate some point sources which are 
blended in the lower resolution SINGS maps. More importantly, the sidelobes of the MIPS 
PSF have been removed around point sources. This is evident in the case of the companion 
and for strong point sources in the disk of M51a. 

HiRes deconvolution is most effective in removing diffraction spikes in the case of isolated 
point sources located in a zero background. However, when the point sources are embedded 
in a high background (like the nucleus of M51a and bright point sources in the spiral arms), 
point sources show a residual ringing at a very low level of 0.2% to 0.5%, although the 
diffraction spikes are removed. This level is significantly higher than the nominal diffraction 
residual around isolated bright point sources. For comparison the ring around the nucleus 
of M51b is weak at the 0.2% level. Therefore the nucleus of M51a and the region around the 
nucleus of M51b will be masked in the correlation analysis applied in the following sections. 
We verified that the flux has been conserved during the deconvolution by comparing the total 
flux of the initial and deconvolved maps, as well as the flux within individual regions. For 
this second test, we used 13" diameter apertures centered on bright H ll regions as defined by 



Kennicutt et al. (2007). The total flux and the fluxes within the apertures centered on the 
bright Hll regions are conserved, with a difference of about 0.7% between the maps and the 
rms noise level of the final HiRes maps are about 0.0067 MJy/sr at 24 fim and 0.0069 MJy/sr 
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at 8 /iin. 



2.3. Ha data 

The Ha image used in our study is the Ha emission-hne map with the underlying stellar 



continuum subtracted as presented by Calzetti et al. (2005). An average constant correction 



factor has also been applied across the image to correct for the contribution of the two 



[Nil] emission lines (Calzetti et al. 2005). It was obtained on 2001 March 28 at the 2.1m 
KPNO telescope as part of the SINGS Legacy Survey. The measured FWHM of the PSF 
is 1'.'9 comparable to the resolution of our radio continuum data. In addition to standard 
reduction procedures, correction for the vignetting effect in the southern region of the FOV 
was apphed to these data, by comparing flux measurement of stars in the vignetted side of 
the mosaic with the same stars in the clean side. Final flux calibration was carried out using 
standard star observations. The final map, without any further background subtraction 
applied, is shown in Fig. [2j The median value of dust extinction in Ha is Ay ~ 2.8 in Hll 



regions ( 


Calzetti et al. 2005 Scoville et al. 


2001 


) and Av ~ 2 


2001 


) and it increases towards the center ( 


Calzetti et al. 


2005) 



2 in diffuse gas (Scoville et al. 



optical light of young stars and re-emit in IR. Therefore, dust extinction hampers greatly 
any SFR estimation from optical and UV wavelengths. 



3. M51's radio continuum as seen by VLA at 100 pc resolution 

3.1. Radio morphology 

Fig. |2] presents our final short-spacing corrected radio continuum maps. These high 
resolution radio images of M51 show unprecedented details. Many structures are revealed 
from point sources in the spiral arms to faint extended filaments in the inter-arm and outer 
regions. While the general distribution of the 20 cm continuum, dominated by synchrotron 
emission, is fairly similar to the morphology seen in the 3.6 cm and 6 cm emission (a mixture 
of synchrotron and thermal radio emission), the contrast of individual structures is changing 
from one frequency to another. 

The spiral arms are bright and well defined at all wavelengths, yet faint extended emis- 
sion is hardly detected at 3.6 cm. Within the arms, compact bright sources are present at 
all three frequencies. They mostly correspond to Hll regions and SNRs. The inter-arm 
regions present filamentary structures as well as fainter extended smooth emission at 20 cm 
and 6 cm. Those structures are not seen at 3.6 cm, due to the higher noise of the 3.6 cm 
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data and its intrinsically higher faintness as evidenced by its spectral index [a ~ —0.9, with 



Sy oc z/", Fletcher et al. 2010). Distinct 20 cm radio continuum emission is present between 



the northern arm of M51a and its companion. This region is less bright in the 6 cm map 
and is hardly evident in the 3.6 cm map. It is consistent with synchrotron emission that is 
enhanced by the ongoing interaction of NGC5195 with M51a (see Schinnerer et al. in prep, 
for a detailed analysis of this region). In particular an interesting feature in this region is the 



arc shaped structure south of M51b, also observed in Ha emission (jGreenawalt et al. 1998 
and Fig. |2 



North of NGC5195, extended faint emission is present at 20 cm and 6 cm. In 
the 20 cm and 6 cm images, one background radio galaxy with strong radio lobes is observed 
east (RA = 13''30"'16."; 6 = +47°10'22") of M51. This object is brighter at 20 cm than 6 cm, 
which implies a steep spectrum source. It is not seen in the 3.6 cm map since it lies outside 
the mosaic. 

M51a hosts a type 2 active nucleus. The inner kiloparcsec region observed at 20 cm 
is shown in Fig. pi The nucleus is unresolved and located at RA = 13'*29'"52.^7; 5 = 



+47°11'42'.'5 (J2000), in good agreement with the position reported by Crane & van der Hulst 



(1992) at 6 cm. In addition to the radio source, two main nuclear features are observed: a 



ring-like structure 5" north of the nucleus, thought to be an outflowing bubble, related to 



the nuclear activity (Maddox et al. 2007) and bright compact emission south of the nucleus 



corresponding to the edge of the radio jet, identified by Crane Sz van der Hulst] ( |1992D at 
6 cm. The jet itself, 3" long, is hardly resolved at our resolution (1.4" at 20 cm) as seen in 
Fig.m 



3.2. Radio fluxes and scale lengths 



Using the task IRING we calculated the total integrated flux in the primary beam 
corrected maps. We used rings centered on the AGN location and assumed a position angle 
of 170° and a disk inclination of 20°. We found total fluxes of 0.3±0.2 Jy at 3.6 cm, 0.6±0.2 Jy 
at 6 cm and 1.4±0.1 Jy at 20 cm (including the companion M51b, see Table. [3]). These fluxes 
are in good agreement with previous single dish and interferometric observation. Indeed at 



3.6 cm Klein & Emerson (1981) derive a flux of 0.306 ± 0.026 Jy, at 6 cm Israel & van der 



cm 


Israel & van der 


(1977 


) give a flux of 



HulstI ( |1983| ) report a flux of 0.525 ± 0.033 Jy and at 20 cm Segalovitz (|1977|) give a flux of 
1.5 ± 0.1 Jy. The global spectral index is a = 

Si, oc u"), in good agreement with the results of [Fletcher et al.| (|2010[). 



-0.9 (where the radio flux is expressed as 



Within each tilted ring, we also derived the azimuthal average radio flux. Fig |4] presents 
the radial proflles at 3.6 cm, 6 cm and 20 cm, from 2'.'5 ~',100 pc to 400" ~ 16 kpc, encompass- 
ing the entire extent of M51. The azimuthal averaged radial proflles are similar at all three 
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wavelengths. They show several peaks at the same location: 25", 67", 120" and 147" (Ikpc, 
2.7 kpc, 4.8 kpc and 5.9 kpc). Those radii correspond to the inner and outer spiral arms. 
One difference between our three maps is that the profile at 20 cm and 3.6 cm decrease for 
radii R> 260" ~ 10.5 kpc and R> 220" ~ 9 kpc, respectively, while at 6 cm the intensity 
remains constant at large radii. The radio profiles at 20 cm and 6 cm present a peak at about 
250" (10 kpc), which is not seen at 3.6cm. This peak corresponds to NGC5195. It does not 
show up in the 3.6 cm radial profile due to increase of noise at these large radii in the 3.6 cm 
map. 

We fitted the radial profiles for radii between 50" and 280" (Fig. |4]) with exponentials 
of the form Iq x exp {—ji), where h is the exponential radial scale length of the profile. 
The corresponding scale lenghts are 3.8 kpc ±0.2 kpc at 3.6 cm, 5.6 kpc ± 0.4 kpc at 6 cm 
and 4.5 kpc ±0.5 kpc at 20 cm (see Table. |3]). Our values are lower than the scale length of 
the molecular gas disk (6.5 kpc; [Hitschfeld et al. 2009) but the scale length at 6 cm is in 



good agreement with the dust scale length (5.4 kpc; Meijerink et al. 2005). At 3.6cm the 
smaller scale length corresponds to the thermal emission, concentrated close to H ll regions. 
We find a smaller scale length at 20 cm than at 6 cm, suggesting that we may miss a very 
extended smooth component, consisting of old CR electrons, emitting synchrotron emission. 
Although our VLA data are theoretically sensitive to structures up to 15' (about 35 kpc), 
single dish observations seem to be needed to correctly map the largest scale disk at 20 cm. 
For comparison, NGC6946 presents a scale length of 4.0 kpc at 20 cm (non-thermal emission) 



and 4.6 kpc at 6 cm (Walsh et al. 2002) while in M33, the scale lengths are much smaller: 



2.7kpc at 20 cm, 2.3kpc at 6 cm and 1.3kpc at 3.6 cm ( Tabatabaei et al.|2007b ). The missing 
flux at 20 cm does not pose a problem for our wavelet analysis, as it only affects spatial scales 
larger than ~400". We verifled this by comparing the wavelet analysis of the 6 cm images 
with and without short spacing correction. 



4. Structural Analysis of the Continuum Morphology 



We want to investigate and quantify how well the structures in M51a at different spatial 
scales (spiral arms, fllaments in the inter-arm regions, clumps in the spiral arms...) correlate 
between the different continuum emission maps, in particular the radio and IR maps. In order 
to study the correlation of structural scales at different wavelengths, the wavelet analysis is 



a powerful tool (see e.g. Frick et al. 2001 Hughes et al.| 2006b Tabatabaei et al. 2007a), as 



it allows for a comparison of compact (e.g. Hll regions), intermediate (e.g. spiral arms) and 
large structures. The mathematical formalism of the wavelet analysis (deflnitions, energy 
conservation and image reconstruction) is described in detail in Appendix |Aj In this section. 
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we present the wavelet analysis of our multiwavelength data of M51. 



4.1. Wavelet Analysis of M51 

To apply the wavelet analysis to our data, we first convolve our radio and IR maps to 
a common resolution of 2! A (which is the resolution of our HiRes 24 /xm and radio 3.6 cm 
images). The maps were regridded to the same geometry with a spatial sampling of (T.'Q 
and centered at RA = 13''29™52."7; Dec = +47°11'42'.'5 (J2000), the location of the nu- 
cleus of M51a. Finally we cut the images to a common field of view of 12'xl2' (799x799 
pixels). We subtracted five bright radio sources as they dominate the wavelet spectra at 
small spatial scales and therefore mask the contribution of fainter features. These sources 
are the two galactic nuclei (including the jet and bubble related to the AGN of M51a and 
the surrounding of the nucleus of M51b) and three luminous radio sources at the coordi- 
nates RA = 13'^29'"30."5; Dec = +47°12'50'.'5, RA = 13^^29^51. "6; Dec = +47°12'08" and 
RA = 13''30'"05"; Dec = +47°10'35'.'7 (J2000). The latter three sources correspond to the 



sources 1, 37 and 104 in the catalog of Maddox et al. (2007). As they exhibit steep radio 



spectral indices and no optical or X-ray counterparts have been detected, they are very likely 



background radio galaxies (Maddox et al. 2007) 



In order to avoid any boundary conditions while scanning the dynamical range of spatial 
scales with the wavelet decomposition, the final maps are 3 times larger than the original 
maps (2397x2397 pixels). We developed a new tool to calculate the wavelet transform using 
MatLab and the YAWTb toolbox (Yet Another Wavelet ToolboxQ, to be able to handle 
such large maps. Within this toolbox, the continuous 2D wavelet transform is implemented 
in terms of Fourier transforms, using the FFT algorithm within MatLab: 



W{a,x) 



f{k)^*{ak)e'''"'-''dk 



— oo J — oo 



where / is the 2D function to be analysed and ijj is the wavelet function. / and ijj are the 
Fourier transforms of / and ip, respectively. 

The choice of the analysing wavelet depends on the data and the goals of the analysis. 



Here we choose the Pet Hat function introduced by Frick et al. (2001) which allows a better 



separation of scales than most commonly used wavelets (e.g. Morlet, Haar or Mexican Hat 
wavelets). The Pet Hat function is defined in the Fourier space as: 



^http://rhea. tele. ucl.ac.be/yawtb/ 
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where k = \k\. 



While the wavelet analysis is robust to noise, compared to Fourier analysis (Frick et al. 



2001 Tabatabaei et al. 2007a), it is not obvious how exactly noise affects the smallest scales 
probed at low signal-to-noise ratio (S/N). In our radio data, the noise at the map edges is 
increased due to the primary beam correction of the interferometric observations. Since we 
have added single dish data to our 6 cm and 3.6 cm VLA maps for short spacing correction, 
blanking pixels of values lower than, for example, 2arms removes also a lot of this intrinsically 
faint, but real, emission. This is not the case for the 20 cm and the IR and optical maps, 
but for consistency we applied the wavelet analysis to our data without blanking pixels with 
low S/N. 

We tested the influence of noise by gradually adding Gaussian noise to our 20 cm map 
(which has the best signal-to-noise ratio of the radio data) from llyuJy/beam (noise level 
of the 20 cm map) to 25/iJy/beam (noise level of the 3.6 cm map) in steps of 1 /xJy/beam. 
The description of the noise effect on the wavelet analysis can be found in Appendix |Bj 
The conclusion of our test is that noise affects the wavelet analysis on small spatial scales: 
a < 18" (a < 10") at the noise level of the 3.6 cm map (the 6 cm map, rms ~ 16 //Jy/beam). 
Therefore, we consider the wavelet spectra of 3.6 cm and its correlation with the other wave- 
lengths only on scales larger than 18", while this limit is set to 10" for the spectrum and 
cross-correlation of the 6 cm data. 

For the final analysis, we decompose the radio, IR and Ha maps into 50 spatial scales 
between 5" (twice the resolution, except for the radio maps at 3.6 cm and 6 cm, see above) 
and 650" (~ 26 kpc, maximum extent of M51a) to encompass all structures within the galaxy 
with enough spatial sampling to compare the morphology at different wavelengths. 



4.2. Dominant spatial scales as revealed by wavelets 

In this section we describe the results of the wavelet analysis of M51a's radio, IR and 
optical maps. Table. |4] summarizes the interesting scales found in the wavelet spectra (labeled 
with a 'S', second and third column) or cross-correlation between two images (labeled with a 
'X', fourth and fifth column). By construction, the steps between the spatial scales are about 
10% of a: 6a f^, increasing toward larger scales. Moreover as it is described in details in 
the following, there are shifts of the interesting scales between the wavelengths. Therefore, 
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for a better clarity of Table. [4], the errors on a listed in the first column are the maximum 
between the steps and the spatial shift of the local extrema discussed below between the 
different emissions. 



4.2.1. Wavelet spectra 

The wavelet spectrum, M(a), represents the distribution of the emitting power as a 
function of the spatial scale a: 

/•oo /"OO 

M(a) = / / \W{a,x)\^dx 



where (a, x) is the wavelet decomposition map at the spatial scales a (see Frick et al. 



2001, and Section A. 2 in this paper for a detailed mathematical definition). Towards larger 
scales, a smooth increase in M(a) is expected if most of the emission is coming from diffuse 
structures, and a decrease if strong point-like or compact structures are present. Peaks in 
the distribution of M(a) indicate dominant scales in the morphology. 

The wavelet spectra derived for various continuum and line maps of M51 are presented 
in Fig. [5] For a better clarity of this figure, the spectra have been scaled by powers of 
10, avoiding any confusion. None of the spectra increases smoothly towards larger scales 
indicating the absence of a dominant diffuse disk. The only general increase is visible at 
20 cm, while the wavelet spectrum at 24 /xm decreases and the spectra of the other frequencies 
remain roughly flat. These general trends are accompanied by some fluctuations, which 
indicates that prominent scales, and thus distinct morphological structures, exist. Local 
maxima and minima in the wavelet spectra of our different maps are listed in Table. |4| in 
the second and third columns respectively. 

All three wavelet spectra of the radio data exhibit similar substructures. The 3.6 cm 
(6 cm) spectrum starts at 18" (10") as the noise dominates at smaller scales (see Appendix[B|). 
The fluctuations in the three radio spectra can be interpreted as several local extrema, 
roughly located at the same spatial scales. All three radio wavelet spectra show a small 
bump at a ~ 20" (0.8 kpc, corresponding to complexes of star forming regions and structures 
within the spiral arms). The second local maximum can be attributed to the width of the 
spiral arms seen in radio and a third maximum appears on a larger scale where the resolved 
structures include the inner disk of M51a, NGC5195 and the outer spiral arms. Interestingly, 
these two maxima appear at slightly larger scales at 20 cm (respectively at a ~ 55" ~ 2.2 kpc 
and 150" ~ 6 kpc) than for the other two radio wavelengths (at about 45" ~ 1.8 kpc and 
~ 135" ~ 5.5 kpc), indicating wider spiral arms and a larger disk at 20cm likely due to the 
propagation of cosmic ray electrons. 
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The most dominant scale in our data is found at a = 10" (0.4 kpc) where the 24 /xm 
emission is considerably strong. No dominant scale is seen in the 8 /im (PAH) or the Ha 
emission. The latter could indicate a strong effect of extinction in the Ha emission by dust 
flattening the spectrum on small scales, (since we would expect a small-scale dominance in 
the Ha emission, similar to the one found for the 24 /xm emission). Small scales (a <20") 
contribute more to the emission at the IR and optical wavelengths than in the radio contin- 
uum. Although the IR emission at 8 /im and 24 /im is mostly related to warm dust, their 
spectra are strikingly different. While the 24 /im emission is strongest at 10", the PAH 8 fim 
spectrum is approximately flat over our spatial scales range. This reflects the fact that the 
24 fim emission is more strongly peaked on star forming regions than the PAH 8 fim. The 
8 fim spectrum reaches a very weak local maximum at about 16" ~ 640 pc. This weak maxi- 
mum could trace the PAHs illuminated by young and massive stars. However PAHs emission 
arising from PDRs surrounding the H II regions ( Helou et al.|2004 Murphy et al.|2006 Bendo 



et al. 2006 ; Gordon et al. 2008 ) is not predominant and the contribution of PAHs heated 



by the interstellar radiation field on larger scales is significant, as shown in previous studies 
(e.g. Bendo et al. 2008). Moreover, dust emission may be important in regions of intense ra- 
diation field at 8 /xm ( Draine fc Li||2007 ) which could also explain the significatn contribution 
of the large scales in the wavelet spectrum of this wavelength. The 8 fim spectrum shows 
a local maximum at a ~ 45" (~ 1.8 kpc) similar to the radio 3.6 cm and 6 cm emission 
spectra (the wavelet spectrum at 20 cm reaches a maximum at 55" ). The Ha and 24 /im 
wavelet spectra present no minimum at this scale, and decrease smoothly out to a ~ 200" 
(~ 8 kpc). This indicates that while structures corresponding to the spiral arms are clearly 
visible at 8 fim they do not dominate the morphology of the Ha and 24 fim emission. Both 
IR and optical wavelet spectra have a weak maximum around scales of 120-130" (~ 4.5 kpc) 
corresponding to structures encompassing M51a inner disk, the outer spiral arms and the 
companion NGC5195. Finally a strong minimum at 200-220"(~ 8.4 kpc) is present in both 
IR and optical wavelet spectra at the same location than the minima in the radio wavelet 
spectra. This is moreover the only obvious feature in the Ha spectrum. This strong local 
minimum in all wavelet spectra may reflect the good contrast between the disks of the two 
interacting galaxies on this spatial scales (~ 8.4 kpc). This minimum is also weaker in the 
wavelet spectra of the radio emission, revealing a diffuse radio structure encompassing both 
M51a and M51b and interestingly the 3.6cm wavelet spectrum reaches this local minimum 
at larger scales (~ 300" ~ 12 kpc). The locations of the local extrema in the optical and 
IR wavelet spectra are similar to those of the extrema present in the radio spectra, albeit 
at systematically shorter spatial scales: the sizes of the corresponding structures are 10% 
smaller than those in the radio maps. This shift between the sizes of optical/IR and radio 
structures is likely reflecting the diffusion of the CR electrons into the ISM which makes the 
radio structures appear wider and larger than the IR and optical counterparts. 
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Finally, at all wavelengths the wavelet spectra reach a maximum around 600" (~ 24 kpc) 
and remain almost constant at this maximum for larger scales. On these scales, the wavelet 
probes the extended galactic disk of M51a. This maximum is reached at smaller scales in 
the IR and Ha than in the radio maps, the difference is about 5 kpc. This shows a larger 
extent of the disk in radio than in IR and Ha, again indicating propagation of CR electrons 



producing the radio emission. This finding is in good agreement with Murphy et al. (2006 



2008 ) who showed that radio images can be interpreted as a smeared version of the IR maps. 



4-2.2. Wavelet cross- correlations 

The wavelet spectrum allows us to identify the dominant energy structures in our data 
and to qualitatively compare them. The wavelet cross-correlation is a useful method to 
compare different images as a function of spatial scales as shown in several studies (e.g. 
Frick et al.|2001 Hughes et al.|2006b ; Tabatabaei et al.|2007a ). Indeed, while standard cross- 



correlation analysis, such as a pixel-to-pixel correlation, can be dominated by bright extended 
regions or large-scale structures, the wavelet cross-correlation r^i allows for the analysis of a 
scale-dependent correlation between two images. It is also very sensitive to any differences 
between two images even on spatial scales where the wavelet spectra appear similar. Such 
differences could be, e.g. an azimuthal off-set between spiral arms or a systematic shift of 



structures of similar size. The definitions and equations used here are detailed in Sect A. 3 



Following Frick et al. (2001 ) we consider the correlation between two images to be significant 
for values higher than r^ > 0.75. Fourth and fifth column of Table. |4]list the local maxima 
and minima in the wavelet cross-correlation of our different maps. 

There is significant correlation between the radio images at all scales larger than ~ 20" 
(Fig. [6]). At all scales, the best correlation occurs between the 20 cm and 6 cm maps while 
the weakest is between 20 cm and 3.6 cm particularly for scales a < 140". Since the observed 
radio continuum is a combination of thermal (free-free) and non-thermal (synchrotron) emis- 
sion, the contributions from thermal and non-thermal emission change as a function of wave- 
length as well as location in the galactic disk, the fraction of non-thermal emission being 
more and more important at long wavelength as well as on large structures. This can ex- 
plain the differences between the wavelet correlations shown in Fig. |6j Two local maxima 
and two local minima are present in these cross-correlations at similar spatial scales and are 
especially pronounced in the cross-correlations involving the 3.6 cm data. The locations of 
the first local maximum and minimum correspond well to those seen in the wavelet spectra 
of the radio maps and corresponds to the width of the spiral arms (about 2 kpc) and the 
size of the structures formed by the combination of the spiral arms and the interarm region 
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(~3.5kpc) . The local minima could be attributed to different slopes and displacements in 
the extrema of the corresponding wavelet spectra (Fig. [s] and Table |4]) again caused by more 
diffuse morphology of the longer wavelength radio emission (due to a higher contribution 
from synchrotron emission). 

Now we study the wavelet correlation between the radio continuum at 20 cm, 6 cm and 
3.6 cm and the IR/optical emission (Figures [7] and [sj respectively). First we notice that 
all the correlation coefficients decrease globally towards smaller scales and this decrease is 
particularly steep from a = 40" downwards, which corresponds to the width of the spiral 
arms. The correlation of the IR emission at 24 /im is higher with the radio continuum at 
3.6 cm than at 6 cm and 20 cm on spatial scales smaller than 100" (Fig. [?]). This is expected, 
as the radio continuum at 3.6 cm contains less synchrotron contribution than at 6 cm or 
20 cm, the high correlation between 3.6 cm and IR data reflects the common origin of the 
thermal radio and IR emission, both being powered by ionizing young stars in the Hll 
regions. Local minima and maxima are present in the three correlations between radio and 
IR emission, appearing at similar spatial scales as the local extrema in the cross-correlation 
between the radio wavelengths themselves. Finally, the correlation between 24 /im and 8 /im 
PAH emission is the strongest correlation involving IR or optical emissions, at all spatial 
scales. It reaches a weak local minimum at a = 245" (~ 7.3 kpc) which is the location of 
the local maximum of the cross-correlation between 24 /im and 3.6 cm. 

The correlations including Ha emission decrease more rapidly than the corresponding 
correlations with IR emission (Fig. |8|. They all present a local maximum at about 135" 
(~5.7kpc) and minimum at 200" (~8.2kpc, this local minimum is reached at slightly larger 
scale for the correlation with 3.6 cm). The wavelet cross-correlations between Ha and the 
IR emissions at 24 /im and 8 /xm are very similar. As in the case of the 24 /im emission (see 
Fig. [T]), the wavelet cross-correlations between Ha and the radio emissions are lower at 20 cm 
than at 6 cm and 3.6 cm at all scales a < 200". 



5. The local radio-IR correlation within M51 

5.1. Characteristic scales 

Here, we discuss important spatial scales in M51a galactic disk, as revealed by the 
wavelet spectra and cross-correlations, and their physical implications. Fig. [9] shows the 
original maps smoothed to 2'.'4 resolution at 20 cm, 6 cm, 3.6 cm, 24 /im, 8 /xm and Ha from 
top to bottom, and the corresponding wavelet decomposition maps at 4 prominent spatial 
scales: 10" (450pc), 45" (1.8kpc), 83" (3.3kpc) and 200" (8kpc). The wavelet decomposed 
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maps show the structures which contribute at the different scales. The wavelet maps at 10" 
(Fig. [9] second column) show small structures which are dominating the wavelet spectrum 
at 24 /im (see Fig. [s]). Those features also contribute largely to the 8//m and Ha emission, 
although their spectra do not show a pronounced peak at this scale. They are formed of 
compact structures (especially in Ha) corresponding to H ll regions and filament-like features 
(especially at 8fim). 

At 45" (~1.8kpc) we distinguish elongated structures defining the spiral arms (Fig. [9| 
third column). This spatial scale is of the order of the typical width of the spiral arms of 
M51a. At this spatial scale the cross-correlations between the radio/IR/optical data reach a 



local maximum (Sect. 4.2.2). It is interesting to note that some filamentary structures in the 
inter-arm regions contribute to the energy at this spatial scale, the length of these filaments 
being of the same order of the width of the spiral arms. 

Finally, at 200"(8kpc) , the wavelet decompositions show the diffuse emission from the 
galactic disk and the region between M51 and its companion (north-east). This interaction 
zone dominates at long radio wavelengths (20 cm and 6 cm), it is weaker at 3.6 cm and 
disappears in the IR and optical maps. This spectral trend indicates that the emission in 
this region is mainly due to non-thermal radiation, such as synchrotron emission. The spatial 
scale corresponds also to a local minimum in the radio correlations between Ha emission and 
the radio and IR emissions, likely due to this extra radio component in the interaction region 
(Fig.§. 

As shown in the previous sections, wavelet analysis is a powerful tool to determine 
the dominant morphological structures in our mult i- wavelengths images. In the case of 
M51a, the scale at a = 83" ~ 3.5 kpc is particularly interesting. It corresponds to a local 
minimum in the wavelet spectrum of all radio and IR images (Fig. |5|, which appears deeper 
at low than high radio frequency, as well as in the correlation between the radio and IR 



maps (see Sect. 4.2.2). This characteristic is also seen in NGC6946 (Frick et al. 2001) and 



M33 (Tabatabaei et al. 2007a) at roughly the same spatial scale, 3.5 kpc. We interpret this 
minimum due to the best contrast between the arm/interarm region. Indeed at this scale, 
the positive part of the wavelet function filters the structures formed by the combination of 
the bright spiral arms and the "darker" interarm regions. The resulting wavelet coefficients 
contain then information from these two regions and on the spatial scale where the contrast 
between spiral arms and interam is optimal, the wavelet coefficients are then lower than at 
smaller or larger scales, for which the contrast is less good and then more energy is contained 
into the positive part of the wavelet function. This leads to a local minimum in the wavelet 
spectrum at the corresponding spatial scale. Thus the wavelet transforms at the scale of 
83" of M51a's radio and IR images trace well the spiral arm structure (see Fig. lol fourth 
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column). We will use this result to investigate the local radio-IR correlation in M51a in 
the following two sections. First we discuss the wavelet cross-correlation between the 20 cm 
and 24 /xm maps, then we use the wavelet transform at 83" to define different regions in 
MSla's galactic disk (spiral arms; interarm region; outer and inner regions) and discuss the 
correlation between 20 cm and 24 /im as a function of the environment. 



5.2. Spiral arms at 20 cm and 24 /xm 

The local radio-IR correlation is formed on one hand by the strong relation between the 
thermal radio (Bremstrahlung) and the warm dust emission present in H ll regions, and by 



the, albeit weaker, correlation between the synchrotron emission and the cool dust (Hoernes 



et al. 1998). Here we discuss the correlation between the 20 cm radio continuum, which is 



the most sensitive tracer for synchrotron radiation in our radio data (Condon 1992) and 
the 24 /im emission within M51a. We use the 24/im emission to investigate the radio-IR 
correlation since it has the same spatial resolution as our radio maps, although it is mostly 
related to warm dust. Synchrotron and IR continuum exhibit many complex structures over 
all spatial scales as shown in Figures[T] and |2} Moreover they do not always arise from the 
same location, and in some regions, non-thermal radio emission is not at all associated with 
IR emission, as in the interaction region between M51a and NGC5195. CR electrons are 
usually produced and accelerated in SNRs and radiate synchrotron emission while propa- 
gating in the magnetized interstellar medium. In the regions where radio emission is not 
spatially associated to IR, synchrotron emission from CR electrons may be enhanced (e.g. 
by strong magnetic field) or produced by a different physical process (such as galactic-scale 
shocks) rather that SF induced shocks. 

The wavelet spectra reveal that the predominant structures at 20 cm and 24 fim are 



different (see Sec. 4.2.1). While small scale structures dominate the IR emission, they carry 
only little energy in the radio continuum maps (Fig. |5]). Those structures correspond to 
bright knots in the spiral arms, associated with Hll regions and stellar clusters (Fig. [9]). The 
spiral arms at 20 cm are wider than at the other wavelengths (~ 60" compared to about 40" 
in 3.6 cm and IR maps as traced by the maximum in the wavelet spectra, see Sect. |4] and 
Table |4]), reflecting the diffuse synchrotron emission extending into the interarm regions. On 
larger scales (83" ~ 3.5kpc, see Fig. [?]) the correlation between 20 cm and 24 /im emission 
reaches a local minimum. The two wavelet transform maps at this spatial scales show similar 
structures, tracing the spiral arms and interarm region (Fig|9| fourth column). However, an 
obvious difference between these two maps is that the spiral arm structure revealed by the 
wavelet transform of the 20 cm emission at this spatial scale is more continuous than the 
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structures tracing the arms in the 24 /im wavelet transform. This could be explained by the 
diffusion of the CR electrons from their acceleration site into the ISM. 

The spatial scale on which the correlation between the NT-radio continuum and the 
IR emission should break is generally assumed to be the diffusion length of CR electrons. 
Following our wavelet analysis, the correlation between 20 cm and 24 fim starts to decrease 
rapidly from scales of 45" (1.8 kpc) downwards. Between 30" and 20" (1.2-0.8 kpc) the 
correlation coefficients get stabilized around a value of 0.54 and then decrease further from 
0.5 - 0.8 kpc towards smaller spatial scales, in good agreement with CR electrons' diffusion 



lengths of about 0.5 kpc in M51 estimated by Murphy et al. (2006, 2008). This shows that 



the radio-IR correlation is not smooth and exhibiting local maxima and minima at different 
spatial scales corresponding to particular morphological structures in the galaxy. This result 
reflects the complexity of the local radio-IR correlation and hints that different processes 
lie at the origin of this relation and that while one process may dominate at a particular 
spatial scales, none is dominating at all scales. In particular, the fraction of non-thermal 
and thermal emission can change with spatial scales, while diffuse IR emission form older 
stellar population can contaminate the emission at 24 /im considered here. These different 
components of the radio and IR relation can explain the observed variations in the radio/IR 
relation (traced here by the correlation between the 20cm and 24 /xm emission). This will 
be developed in details in a further study, where in particular careful decomposition of the 
radio continuum into thermal and non-thermal components will be done (Dumas et al., in 
preparation). The decrease of the correlation coefficients towards smaller spatial scales can 
be interpreted as a break of correlation between the radio and IR emission on those scales. In 
particular it is interesting to note that the correlation between the 20cm and 24/xm emission 
does not hold on typical scales of star forming complexes. The break of correlation on these 
small scales is expected given the timescales difference between the IR and non-thermal radio 
emission and the diffusion length of cosmic rays in the ISM. 



5.3. Observed environmental effects 

In order to investigate the radio-IR correlation in the different environments of M51a, we 
identified different regions within the galactic disk. The spiral arms and interarm region are 
defined by the wavelet transform at the spatial scale a = 83", as this spatial scale provides 



the best contrast between the arm/interarm region (see Sec. 4.2.1). Positive values of the 
wavelet filtered 20 cm image at 83" define the spiral arms, while negative values mask the 
interarm region. The region with R< 1.75 kpc is the inner region and the outer disk region 
starts at the outer edges of the external spiral arms, at roughly R> 6.5 kpc. We excluded 
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the companion NGC5195 from the analysis. The different regions are filled with apertures 



of 6" diameter. Figure 10 shows the different regions overlaid on the 20 cm map. The 20 cm 



and 24 /im flux densities within these apertures are calculated with the task 'aphot' in IRAF. 



Fig. 11 presents the flux densities at 24 /xm and 20cm in the spiral arms (top left panel), 
interarm (top right panel), outer disk (bottom left panel) and inner region (bottom right 
panel) . 



We fitted the relation between the 24 //m and 20 cm emission with the function 



24 fim 



b X {S20cmY, with F24^m and S20cm being the flux density at 24 /xm and at 20cm (in Jy), 
respectively. The fitted values and standard deviations of the scaling exponent c and the 
coefficient b, are listed in Table. |5} Prior to this fit, some points in the different regions were 



excluded. Those points are marked in Fig. 12 They correspond to AGN related emission (in 



particular the radio jet) and bright Hll regions in the inner region (white circles in Fig. 12), 
to bright H ll regions in the spiral arms and interarm region (respectively, red and blue circles 



in Fig. 12 ) and to one point at very low 24 /xm flux density in the outer region (green circle 



in Fig. 12). Those outliers are also highlighted as green squares on Fig. 11 The results of 



the fit in each region are given in Table |5] and plotted as red lines in Fig. 11 

The relation between F2A^m and 5*20 cm is linear in the spiral arms and globally over the 

galaxy. In this case, the logarithm of parameter b is equal to the parameter g24 defined by 



Murphy et al. 



(2006) as g24 = log 



24 iim 



{Jy) 



and we find that our fitted parameter is 



in good agreement with the values derived by Murphy et al. ( 2006[ ) (see Table [s]). In the 



interarm region and outer disk, the flux density at 24 /zm is proportional to about the square- 
root of the 20 cm flux density while in the inner regions we observe an opposite behavior: 
F2iij,m oc (5'20cm)^'^- In thcsc regions where a non-linear relation occurs between F2i^m and 
S20cm, 924 is determined by both log(6) and c, hence is not constant and of a limited use as 
a quantitative measure of the radio-IR correlation. 



The large differences in the exponent < c > of the infrared-radio correlation (Table |5 
demonstrate that the physical processes involved are different in the arm, interarm and 
inner regions, while the interarm and outer regions show a similar behaviour. As the spiral 
arms dominate the total flux emission (50% at 20cm and 60% at 24 /im), a linear behaviour 
is obtained when looking at the whole galaxy, as observed for the integrated emission of 



most galaxies (Yun et al. 2001 Bell 2003). Our results are especially striking in view of the 



many processes involved: star formation, dust heating, dust opacity, cosmic-ray acceleration. 



cosmic-ray propagation and magnetic fields (see e.g. Bicay & Helou 1990 ; Murphy et al. 2006 ). 
In spite of the complicated interplay of the ISM's different components, three well-defined 
regimes can be distinguished and are discussed in detail below. 
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5.4. Possible causes for environmental dependencies 



Here, we are following the formalism and assumptions used by Niklas & Beck (1997). 



Assuming equipartition between the energy densities of the magnetic field and CR electrons, 
we can couple the radio non-thermal fiux density to the gas volume density. The energy 
density equipartition assumes that the energy densities of cosmic rays (electrons and protons) 
and magnetic field are approximately equal and gives: 



s. 



20 cm 



oc B 



eq 



where Beg is the total equipartition magnetic field strength (Beck &; Krause 2005) and a 



is the synchrotron spectral index. The equipartition magnetic field strength B^q is also 
referred to as the minimum energy magnetic field as it is the minimum value of B to realize 



this energy equipartition. This assumption has been and is still discussed (Niklas & Beck 



1997 Beck & Krause 2005; Thompson et al. 2006). In particular Thompson et al. (2006) 



compare this value of the magnetic filed strength measured in different galaxies to different 
scaling relations between B and the gas surface density and show that the equipartition 
magnetic field strength underestimates the magnetic field in starburst galaxies. M51 however 
follows the relation of normal star forming galaxies for which the magnetic energy density 
of the equipartition magnetic field is comparable to the total hydrostatic pressure of the 
ISM (see Eq. 3 and Fig. 1 of Thompson et al. 2006). In that case we can assume that the 



equipartition magnetic field strength refiect to total magnetic field in M51 and that energy 



density equipartition is a good assumption for this galaxy. In M51 Fletcher et al. (2010) find 
a ~ 0.6 in the spiral arms and inner region, and a ^ 1.2 in the interarm and outer regions. 
In addition, the magnetic field is coupled to the gas: 



B oc pj, 

(3 describes the coupling between the magnetic field and the gas volume density Pg. For 
turbulent field amplification (e.g. by dynamo action), /3 ~ 0.5 is expected (Beck et al. 



1996 Thompson et al. 2006). This value also emerges from this time equipartition between 
magnetic and turbulent energy densities. A similar value is predicted for cloud collapse 



along magnetic field lines (Mouschovias 1976), while isotropic collapse in turbulent fields 
gives /3 ^ 2/3. In case of one-dimensional compression or shear of magnetic fields by 
large-scale fiows of diffuse gas, /3 ~ 1 is expected. 

We therefore obtain a relation between the radio NT emission and the gas volume 
density: 

>->20 cm lX. Pg 
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On the other hand we can also hnk the 24 /im flux density to the gas density. In the case 
of optically thick dust for UV photons, the dust re-processes the totality of the radiation of 
the ionizing stars and thus the IR emission traces the star formation rate (SFR). Therefore 
we assume: -^24^^ -^0 where Jq is the energy density of the dust-heating radiation field 
and Jo oc SFR where SFR is the star-formation rate, hence ^24^™, oc SFR. Moreover the 



Schmidt law (Schmidt 1959) gives a relation between the star formation rate SFR and the 



gas volume density pg. SFR oc p^, with n being the Schmidt law index. Therefore we can 
write: 

F24:fim oc Pg 

and n 



and 



1.4 ± 0.3 (Niklas & Beck 1997). In that case, the correlation between ^24^^ 
5*20 cm with the exponent c yields 

/3 = TTTT^ (3) 



c(3 + a) 



If the dust is not optically thick for UV photons, as, e.g., in the galaxy M31 (Hoernes 



et al. 1998), then the optical depth of the dust r has to been taken into account: 



24 fim 



OC r/o 



Jo is the energy density of the dust-heating radiation field directly proportional to the star 
formation rate, and r oc Pdust, then assuming a constant dust-to-gas ratio we have a relation 
between the 24 pm IR flux and the gas: 

-^24/^01 OC p 



and Eq. [3] becomes 



/5 



n 



c(3 + a) 



(4) 



As shown in Calzetti et al. (2005), in most regions of M51, (e.g. the spiral arms) dust 
attenuation is high enough that the IR emission traces well the SFR and Eq. |3]can be applied. 
However, in region where dust extinction is lower and diffuse IR emission, coming from dust 
heated by the old stellar component (the so-called 'Cirrus' emission) dominates, as at large 
radii or in the interarm region, IR emission is not totally related to recent star formation. 
The IR emission for a fixed SFR will be lower in such regions. And in the case of diffuse 
Cirrus emission we have the following relations: Fdrms oc pdust oc pg , where pdust and pg 
are the volume density of dust and gas respectively (we assume a constant dust-to-gas ratio). 



Then, ^24^^ oc Pg and Eq. [3] becomes 



/9 



1 

a(3 -|- a) 



(5) 
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We have derived three possible relations between the IR and 20cm emission, depending 
on the properties of the ISM. We now discuss the three trends found between the radio and 
IR emission in the spiral arms, interarm region, outer disk and inner region of M51a: 



1. In the spiral arms, we can assume the dust to be optically thick and directly apply 
Therefore, with < c >= 1.01 ± 0.01, a = 0.6 ± 0.02 ( [Fletcher et "aL][2Moj ) 
= 1.4 ± 0.3, we obtain /3 = 0.4 ± 0.1 which is slightly below the range of 



Eq.[3| 

and n 



theoretical predictions, but with an overall good agreement. The expected value under 
the assumption that magnetic density energy is comparable to the turbulent energy 
density of the ISM is /3 = 0.5 (iBeck et al.|p96| [Thompson et al.|[2006|. Fig. 1 



of Thompson et al. (2006) shows that the equipartition magnetic field strength in 



M51 realizes this assumption., and in that case, the Schmidt law index should be 
n = 1.8 ± 0.2. 



2. In the interarm and outer region, the spectral index is higher a 
et al.||2010p. Therefore, with < c >= 0.51 ± 0.01, and n = 



1.2 ± 0.06 (Fletcher 



1.4 ± 0.3, we obtain 

(3 = 0.7 ±0.1 in the case of optically thick dust. Assuming UV-thin dust for the 
interarm and outer regions, we obtain /3 = 1.0 ± 0.1, as expected for compression 
or shear of magnetic fields in flows of diffuse gas. This possibility should be further 
investigated by comparing high-resolution data in radio continuum, H I and CO. 

In addtion, the Cirrus IR emission may provide another interpretation of the non- 
linearity between F24^m and 5*20 cm in the interarm and outer regions. In this region, 
we find that the radio-IR correlation is sub-linear, which implies that for a fixed SFR 
the 24 /xm is lower that predicted (assuming 20cm traces the SFR), as expected for 
regions dominated by Cirrus IR emission. Applying Eq. |5] we then find /3 ~ 0.5 in 
the interarm and outer regions, as in the spiral arms, which is in good agreement with 



theoretical values from turbulent field amplification (Beck et al. 1996). Therefore 



diffuse IR emission emitted by dust heated in large part by older stars seem a good 
explanation for the sub-linear radio-IR correlation found in the interarm and outer 
regions. 

Further high-resolution observations and measurements of dust opacity and separation 
of the cold and warm dust components are needed to distinguish between these two 
possibilities. 

3. For the inner region we found < c >= 1.5 ± 0.1. In this region, the AGN related 
radiation, in particular the radio jet, has been excluded for the fit. With a ^ 0.6 
([Fletcher et"aI]|2010D and n = 1.4 ± 0.3 we obtain /3 = 0.26 ± 0.05, in the case of 



optically thick dust. This value of /3 is well below the range of theoretical predictions. 
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The equipartition field is too weak compared to the gas density. Either the process 
amphfying the magnetic field is less efficient in the inner region, or the equipartition 
assumption yields field strengths which are too low. The former case needs further 
investigation which is beyond the scope of this paper. The latter case is known to 
be valid for strong synchrotron or inverse Compton losses of the cosmic-ray electrons 

. As the magnetic field and the radiation field have large energy 
densities in the inner region of M51a, the radio synchrotron emission is probably sup- 
pressed, hence the radio-derived SFR underestimate the true amount of star formation 
present in the inner 1.5 kpc. 

In summary, we can explain the different exponents of the infrared-radio correlation seen 
in the different galactic environments by a change in the optical dust opacity and the rela- 
tions between the magnetic field strength and the gas density and the presence of a 'Cirrus' 
contribution to the IR emission. The steep exponent detected in the inner region is probably 
the effect of suppression of radio synchrotron emission due to energy losses of cosmic-ray 
electrons. Future high-resolution observations of galactic disks with HERSCHEL and the 
EVLA (and later with the SKA) will provide an excellent tool to study the interactions of 
the ISM components in galaxies. 



(Beck fc Krause1[2005 



6. Conclusion 

In this paper we presented high quality radio continuum data of the nearby galaxy 
system M51 at three frequencies, 1.4 GHz (20 cm), 4.9 GHz (6 cm) and 8.4 GHz (3.6 cm), at 
a resolution between 1.5" and 2.4". Combining these radio data with deconvolved Spitzer 
IR images, we were able to investigate the radio-IR correlation within M51 on spatial scales 
down to about 200 pc. The wavelet analysis of the radio and IR maps of M51 reveals 
that the radio-IR correlation is not uniform accros the galactic disk and exhibits a complex 
behavior due to the different galactic structures present in the galactic disk. In particular, 
the radio-IR correlation is high within the spiral arms while it becomes weaker for scales 
encompassing the spiral arms and interarm regions, indicating that another process unrelated 
to star formation might contribute to the radio emission, such as a local enhancement of the 
synchrotron emission by the magnetic field. 

Using our wavelet analysis to define different environments within the disk of M51 we 
found that the IR-radio correlation is globally linear over the galaxy and within the spiral 
arms, in apertures of 6" (240 pc) diameter, reffecting that both IR and radio emission are good 
SF tracers in regions where dust exctinction is high. For the first time we could quantitatively 
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show that the radio continuum-IR correlation is not hnear in the interarm and outer regions, 
where the 24 /xm emission is proportional to the square-root of the 20 cm radio continuum 
emission, as well as in the central part of M51 where the radio-IR correlation is over- linear. 
This can be well explained by variation of the dust opacity and the predominance of a diffuse 
IR emission coming from the old stellar population in the outer and interam regions, while 
the validity of estimates of magnetic field strengths based on the equipartition assumption 
are reflected in the exponent of the correlation in the inner regions. 
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A. Wavelet analysis: mathematical formalism 

A.l. The 2-D continuous wavelet transform 

The wavelet analysis of an image is a generalisation of the Fourier analysis, where the 
sinusoidal waves are replaced by a family of self-similar functions that depend on scale and 
location. These functions are generated by dilating and translating and rotating a mother 



function which is called the analysing wavelet. Frick et al. (2001) showed that the wavelet 
analysis is well adapted for the investigation of scale-dependent properties in astrophysical 
objects. For this study, we restrict our analysis to the use of isotropic wavelets, we assume 
here that the structures that we will probe with such analysis are axisymmetric, hence with 
no favored orientation. We then do not apply any rotation to the mother wavelet. The 
continuous 2-dimensional wavelet transform can be written as: 

W{a,x) = co{a) H H fix')tP* f ^-1^^ dx' (Al) 



— oo ^ — oo 



where a is the scale size, / is a 2-dimensional function (the image), ip* is the complex 
conjugate of the analysing wavelet, x and x' are coordinate vectors and u{a) is a normalisa- 
tion factor. In the following this factor will be equal to ^ for energy normalisation between 
the different scales, the studied images and the wavelet transforms at different scales being 
then in the same unit. 



A. 2. Wavelet spectra 



Following Frick et al. (2001) we defined the energy of the wavelet decomposition map 
at scale a as: 



M(a) 



oo POO 



\W{a,x)\'^dx 



(A2) 



— oo J — oo 



where W{a, x) is the wavelet decomposition map at the spatial scales a. 

The distribution of this energy with the scale a is the wavelet spectrum. This spectrum 
reveals the dominant scales in the initial image. 
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A. 3. Wavelet cross-correlation 



The cross-correlation is a useful method to compare different images. Frick et al. (2001) 



showed that the pixel-to-pixel cross-correlation is dominated by the large scale structures. 
The information given by the wavelet transformations is used to probe the correlation be- 
tween different wavelengths, for example, on the scales of the spiral arms, the interarm and 



central regions. The wavelet cross-correlation coefficient at scale a is defined following Frick 



et al. (2001): 



r^(a) 



JJWi{a, x)W2{a, x)dx 
[Mi(a)M2(a)]^/2 



(A3) 



where Wi and W2 refer to the wavelet decomposition maps at scale a of two images of 
same size and geometry and Mi and M2 correspond to the wavelet coefficients at this scale 



see Eq. A2), and the uncertainty on this cross-correlation is defined as: 



Ar,„fa) 



v/l-rg,(a) 



(A4) 



with = (L/a)^, L being the linear size of the image. In this analysis we use r2 = 0.75 



as the arbitrary threshold for a significant correlation, following Frick et al. (2001) and 



Tabatabaei et al. (2007a) 



A. 4. Energy conservation 

The wavelet transform is conservative: there is no loss of energy when applying the 
wavelet transform to an image f{x). We indeed have a Plancherel relation for the wavelet 
transform, similar to the one for the Fourier transform, in the case of an axisymmetric 
wavelet: 

||/|P = Jf{xrdx=^jM{a)^-^ (A5) 

where M(a) is the energy of the wavelet transform of f{x) at the spatial scales a and 

Ci/i = f — 7T, — dk a constant, with xjj the wavelet Fourier transform. This conservation 
equation is exact for continuous scales from to infinity. Of course, in practice we need 
to choose a finite numbers of scales on which the wavelet transform is produce. However, 



Equation X5 can be used to quantify the contribution of the structures between 2 spatial 
scales to the total energy. 



A. 5. Synthesis 

As a consequence of the energy conservation, we have an exact reconstruction of the 
image / from its wavelet decomposition: 

/(^) = ^ / W{a, h)i. f ^) rfb^ (A6) 
C ^p J \ a J a 

with W{a,h) being the wavelet transform of f{x), with respect to the wavelet ip. As 
for the energy conservation equation, this reconstruction is mathematically exact but not 
realizable in practice since it requires an infinity number of scales from to +oo. Image 
reconstruction (and compression) with wavelets uses the discrete wavelet transform, for which 
an exact image reconstruction can be achieve with a finite and discreet number of scales. 



B. Wavelet analysis: effect of the noise 



We tested the influence of noise by gradually adding Gaussian noise to our 20 cm map 
(which has the best signal-to- noise ratio of the radio data) from 11 /iJy/beam (noise level of 
the 20 cm map) to 25 /xJy/beam (noise level of the 3.6 cm map) in steps of 1 /iJy/beam. We 
then compute the corresponding wavelet spectra, M(a), for the 20 cm map and each of the 
noisy maps, and the wavelet cross-correlations between the 20 cm and noisy maps (see 
Sect. 4.2 and Appendix |A] for the definitions of the wavelet spectrum and cross-correlation). 



The result of this test is shown in Figure [13} where we plotted only the curves corresponding 
to the 20 cm , 6 cm and 3.6 cm noise levels, for clarity. We define the relative differences 

between the initial and noisy maps as AM(a) = ^ 



rw, noisy (a), where M, 



init 



(M, 



noisy J 



^noisy(,(^) 



and Ar„ 



1 - 



Minitia) 

is the spectrum of the initial 20 cm map (each noisy map) 
is the wavelet cross-correlation between the initial 20 cm map and each noisy 
map. The value 1 in the expression of Ar^ correspond to the maximum in the correlation. 



and r. 



w,noisy 



Fig. 13, (top left) shows that the wavelet spectra of the noisy maps differ at small scales, 
as the noise affects the spectra by adding energy at those scales. On the other hand, the 
cross-correlations between the initial and noisy 20 cm maps show that if the noise increases. 



the cross-correlation decreases rapidly (Fig. 13, top right). The relative differences between 



the spectra and the cross-correlations are larger than 10% at scales smaller than 18" for the 
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noisiest map, corresponding to the noise level of the 3.6 cm map. At the noise level of the 6 
cm map (16//Jy/beam) the relative differences reach 10% for scales smaller than 10". 



33 




IS^SO^lS^ 30^00=^ S9"'45= 2U"'J0" 13"U0"'IU' 30™00= a9"'45= 29'"30= 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 1. — SINGS (left) and HiRes deconvolved (right) images of Spitzer MIPS 24 /xm (top) 
and PAH IRAC 8/xm (bottom) observations of M51. The horizontal lines in the top right 
panel correspond to 10"~0.4kpc, size of complexes of star forming regions, and to 83"~3.3kpc 
corresponding to the scales encompassing the width of the spiral arms and interarm regions. 
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Right Ascension (J2000) Right Ascension (J2000) 

Fig. 2. — Radio continuum images at 20cm (top left), 6cm (top right) and 3.6cm (bottom 
left) and the Ha map (bottom right). The data at 3.6 cm and 6 cm have been short-spacing 
corrected using single dish data from the Effelsberg 100 m telescope and all radio maps have 
been corrected from primary beam attenuation. The resolution of the maps is given in 
Table [l] The horizontal lines correspond to 10"~0.4kpc (bottom right panel), size of com- 
plexes of star forming regions, and to 83"~3.3kpc corresponding to the scales encompassing 
the width of the spiral arms and interarm regions (top left panel) . 
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Right Ascension (J2000) 



Fig. 3. — Radio continuum image at 20 cm of the nuclear region of M51. Contour levels are 
at 20(7, 40(7, 60(7, 100(7 and 200(7, with a — 11 /xJy beam"^. The positions of the nucleus 
and radio jet are labeled. 
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10 I ' ' ' ' ' ' ' ' ^ 

50 100 150 200 250 300 350 400 450 
r (arcsec) 

Fig. 4. — Radial profiles of the radio flux at 3.6 cm (top), 6 cm (middle) and 20 cm (bottom). 
The straight lines correspond to an exponential fit of the profiles. The peaks at 25", 67", 
120" and 147" correspond to the spiral arms (marked by 'SP'), the companion NGC5195 is 
seen as a peak at 250"(marked by 'C'). 
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Fig. 5. — Wavelet spectra of the 20cm (blue circles), 6cm (red squares), 3.6cm (black 
pluses). Ha (pink triangles), 24 fim (green diamonds) and 8 fim (cyan crosses) images at 2.4" 
resolution. The spectra are shown in arbitrary units. The scales of the x-axis are presented 
in arcsec (bottom) and kpc (top). The arrows indicate the positions of some interesting 
scales, discussed in the Sections |4.2.1 and 5A_ 
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Fig. 6. — Wavelet correlation spectra between images at 20 cm and 6cm (red circles), 20cm 
and 3.6 cm (blue triangles) and 6 cm and 3.6 cm (green squares). The error bars have been 



computed following Frick et al. (2001) (see Appendix A.3). The horizontal line corresponds 



to a correlation of = 0.75, assumed to define significant correlation above this line (see 
text for details). The arrows indicate the positions of some interesting scales, discussed in 
Section |4X2] and O 
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Fig. 7. — Wavelet correlation spectra between 20cm and 24 /xm (blue diamonds), 6cm and 
24 yum (green triangles), 3.6 cm and 24 fim (red squares) and 8 /im and 24 /im (black crosses). 



The arrows indicate the positions of some interesting scales, discussed in Section 4.2.2 and 

o 
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a (arcsec) 

Fig. 8. — Wavelet correlation spectra between 20cm and Ha (blue diamonds), 6 cm and Ha 
(green triangles), 3.6cm and Ha (red squares), 24 /xm and Ha (black circles) and 8 fim and 
Ha (cyan crosses). The arrows indicate the positions of some interesting scales, discussed in 
Section liSJ 
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Fig. 9. — Original maps (left) of the 20 cm, 6 cm, 3.6 cm, 24yum, 8 fim and Ha emissions 
of M51 (from top to bottom) and the corresponding wavelet decomposition maps at four 
selected spatial scales. The field of view is 9'Axll'.2. 
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Fig. 10. — 6" diameter apertures plotted on the 20 cm image. The apertures fill the different 
regions studied: the central region (white), the spiral arms (red), the inter-arm region (blue) 
and the outer disk (green). 
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Fig. 11. — 24 /iin flux density plotted against 20 cm flux density computed within 6" diameter 
apertures for four regions defined in the galactic disk (see text): spiral arms (top left), 
interarm region (top right), outer region (bottom left) and inner region (bottom right). The 
red lines correspond to the fit derived here (see Tab. [s]). Outliers excluded from the fits are 
marked in green. 



Fig. 12. — Positions of the outlier points which have been removed prior to the fit between 
24 /xm and 20 cm fluxes in the inner region (white circles), spiral arms (red circles), interarm 
region (blue circles) and outer disk (green circle). 
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Fig. 13. — Left: wavelet spectra M(a) of the 20 cm maps with different noise levels (top) 
and the relative differences between the spectra of the initial 20 cm map and each noisy map 
AM(a) (bottom). Right: wavelet cross-correlations between the initial and noisy 20cm 
data (top) and their relative distance to 1, the perfect correlation, Ar^, (bottom ). The red 
and blue lines in all four panels correspond to the noise level of the 6 cm and 3.6 cm data, 
respectively. The horizontal line in the two bottom plots corresponds to a relative difference 
of 10%. 
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Table 1. Data sets used in this study. 



A 


telescope 


resolution 

["x"] 


rms 


20 cm 


VLA (ABCD) 


1.4x1.3 


11 /iJy/beam 


6 cm 


VLA (BCD) + Effelsberg 


2.0x2.0 


16 /iJy/beam 


3.6 cm 


VLA (CD) + Effelsberg 


2.4x2.4 


25 /iJy/beam 




SST/MIPS HiRes 


2.4x2.3 


0.0067 MJy/sr 


8 ^iui PAH ^ 


SST/IRAC HiRes 


1.3x1.1 


0.0069 MJy/sr 




2.1m KPNO 


1.9x1.9 


1.8 X 10^^''ergs~"^ cm~^ arcsec^^ 



^data from the SINGS DR5 
'^data from Regan et al. (20061 
'^data from Calzetti et al. (2005) 
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Table 3. Radio properties of M51. 



A 


20 cm 


6 cm 


3.6 cm 


Su (Jy) 


1.4 ±0.1 


0.6 ±0.2 


0.3 ±0.1 


h (kpc) 


4.5±0.5 


5.6±0.4 


3.7±0.2 


Jo (/iJy/beam) 


55±9 


30±4 


37±5 



Note. — Values were measured using the task IRING in AIPS. The total flux S^, is 
measured out to a radius of 400". The scale lengths h and origin point Jq come from fits to 
the radial profiles presented in Fig. 4 with exponentials of the form Iq x exp (— 
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Table 5. Results of the fit: -^24^^ = b{S20cmT 



fit parameters 


spiral arms 


interarm 


outer 


inner 


total 


(c) 


1.01 


0.51 


0.48 


1.5 


1.00 




0.01 


0.01 


0.01 


0.1 


0.01 


(log b) 


0.9 


-0.50 


-1.19 


2.3 


0.81 


0"log6 


0.01 


0.01 


0.02 


0.3 


0.03 


(^24) 


1.01 


0.91 


0.63 


1.0 


0.86 


'^(124 


0.2 


0.16 


0.17 


0.1 


0.3 



Note. — Means and dispersions of the fit parameters b and c for the different regions, 
assuming F2ifj.m = b{S20cmY (fluxes derived within 6" apertures). The means and dispersions 
of the 024 parameter from Murphy et al. (2006) are also listed for comparison. 



